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1 .  Introduction 

Although  the  Feynman  path- integral  formulation1-  offers  a  general 
approach  for  treating  quantum-mechanical  systems,  only  several  problems  can 

be  solved  exactly.  Two  of  these  are  the  driven  harmonic  oscillator  with  a 

2 

quadratic  Hamiltonian  and  the  time -dependent  damped  driven  harmonic 

3 

oscillator.  A  number  of  situations  such  as  superconducting  quantum 

4  5 

interference  devices,  quantum  nondemolition  measurments, 

magnetohydrodynamics , ^  etc.,  can  be  described  by  driven  copuled  harmonic 

oscillators.  Introducing  the  Caldirola-Kanai  Hamiltonian,^  one  can  obtain 

the  time -dependent  Schroedinger  equation  for  the  damped  harmonic  oscillator. 

However,  it  has  been  a  matter  of  debate  as  to  whether  or  not  this 

8 

Schroedinger  equation  represents  the  quantum  mechanical  dissipative  system. 

9  10 

Some  workers  claim  affirmation  while  others  object  to  it.  This  problem 

11  12 

has  been  critically  reviewed  by  Greenberger  and  Cervero  and  Villaroel. 

The  purpose  of  this  paper  is  to  derive  the  propagator  for  a  driven 

13 

coupled  harmonic  oscillators  (DCHO)  system  from  our  previous  work  for  both 
coupled  and  coupled  driven  harmonic  oscillators  by  means  of  the  path- 
integral  method.  We  introduce  two  harmonic  oscillators  that  are  coupled 
together  with  another  spring.  We  review  the  classical  case  and  construct 
the  form  of  the  propagator  for  DCHO,  respectively,  in  Secs.  2  and  3. 

Section  4  gives  the  exact  derivation  of  the  propagator  for  the  coupled 


harmonic  oscillators  (CHO),  and  in  Sec.  5  we  evaluate  the  exact  propagator 
for  DCHO  by  using  the  results  obtained  in  Sec.  4.  The  energy  expectation 
values  of  CHO  are  evaluated  in  Sec.  6,  and  finally  we  give  results  and 


discussion  in  Sec.  7. 


ii 


By  - 

Distribution/ 
Availability  Codes 
| Avail  and/or 
Dlst  I  Special 


□  D 


2 .  Classical  case 


In  this  section  we  consider  a  system  of  two  harmonic  oscillators  which 
are  coupled  together  by  means  of  another  spring.  We  assume  that  the  masses 
of  the  oscillators  and  three  spring  constants  are  all  the  same .  Let  the 
forces  f^(t)  and  f2(t)  exerted  on  the  two  oscillators  and  their 
displacements  be  and  Xj.  Then  the  Hamiltonian  for  DCHO  can  be  written  as 

H  -  ^  (Pi+P^  +  “^(x^-x^+x2)  *  f1(t)x1  -  f2(t)x2  ,  (2.1) 

2 

where  u>  -  k/m.  Hamilton's  equations  of  motion  for  Eq.  (2.1)  are 


-  Pj/m  (.2.2) 

x2  -  p2/m  (2.3) 
px  -  mw2(x2-2x1)  +•  f1(t)  (2.4) 
P2  -  mw2(x1*2x2)  +  f2(t)  (2.5) 


Equations  (2.1)-(2.5)  yield  the  Lagrangian, 

L  -  (PjXj+PjXj)  -  H 

-  ^  (x2+x2)  -  mw2(x2-x1x2+x2)  +  f1(t)x1  +  f2(t)x2  ,  (2.6) 

with  the  corresponding  equations  of  motion 


(2.7) 


2 

x^  +  u>  (2x^-x2)  -  f^(t)/m 

x2  +  w2(2x2*x1)  -  f2(t)/m  (2.8) 

The  classical  solutions  of  Eqs .  (2.7)  and  (2.8)  are  given  by 

x^(t)  -  A  sin(a>t)  +  B  cos(wt)  +  C  sin(/3wt)  +  D  cos(/3o>t) 

+  J  dr  f  dv  e^w^2r  v  +  f 2 ( )  1  ,  (2.9) 


and 


x2(t)  -  A  sin(cjt)  +  B  cos(wt)  -  C  sin(/3wt)  -  D  cos(/3wt) 

+  |  dr  J  dv  elu<2r-l/-c)  -  f2(i/))  (2.10) 

3 .  Path  integral  of  driven  coupled  hanonic  oscillators 

In  the  path- integral  formulation,  the  solution  of  the  Schroedinger 
equation  is  given  as  the  path -dependent  integral  equations  with  propagator 
K, 


^(x1,x2,t)  -  J  dx^dx^  K(x1,x2,t;  x^,x£,0)  V(x^,x£,0) 


(3.1) 


which  gives  the  wavefunction  ^(x^,x2>t)  at  time  t  in  terms  of  the  wave 
function  0(x£,x2)  at  time  t  -  0.  The  propagator  in  Eq.  (3.1)  can  be  written 
by  means  of  the  Feynman  path  integral 


5 


K(x^,x2,t;x^,x2,0) 


(x1,x2,t) 

(x|,x£,0) 


Dx(t) 


exp [ ( i/Jrf)  S (xx , x2 , , x£ ; 


t)  (3.2) 


where 


1  N‘L 

Dx(t)  -  lim  j  n 
N-KO  A  j-1 


[dXijdXzj/A2] 


(3.3) 


and  S(x1,x2,xi .xi;  C)  is  the  action  defined  as 
Lagrangian  L(x^ , x2 , x^ , x2 ;  t)  between  t  -  t  and 


the  time  integral  over  the 

t  -  0:1 


5(x1,x2,x£,xj;t)  -  J  dt  L(xl,x2,x1,x2;  t) 


(3.4) 


In  Eq.  (3.3)  A  is  the  normalization  factor  given  by 

A  -  [2wi)rfe/m]  ,  e  -  lim( t/N)  (3.5) 

N-«o 


Substituting  Eq.  (2.6)  into  Eq.  (3.4),  the  action  becomes 


S (x^ , x2 , x^ , x2 ,  t)  ■  S ^ (x^ » x2  , x^. , x2  \  t) 


I0  dr  2(yi(°  +  'y\{T)  '  ^[y^CO  -  y1(r)y2(r)  +  y2 ( r )  )  ) 


(3.6) 


where  Sc  is  the  classical  action  and  y^  is  the  deviation  of  x^(t)  from  its 
classical  path  x  ^  given  as 


6 


7l  -  x.  -  xci  .  (i  -  1,2) 


(3.7) 


Then  the  propagator  [Eq.  (3.2)]  can  be  expressed  as 


K(x1,x2,t;  x^.x^.O)  -  F(t)  e 


isc/X 


(3.8) 


Here,  F(t)  is  the  multiplicative  function  given  in  the  form 


-f 


F(t)  -  Dx(t)  {exp  [ ( im/2)rf) 


[  dt  (y?  +  y, 
J0  L 


2  2 

2«  (y[ 


yiy2  +  y2]> 

(3.9) 


It  is  easy  to  show  that  F(t)  has  the  same  form  for  CHO  and  DCHO.  Therfore, 
the  propagator  depends  only  on  the  classical  action  in  both  cases.  In  Eq. 
(3.9),  change  the  variables  x^  ±  x2  into 


Z1  “ 


n 


x2) 


-  72  (x 


1  +  x2> 


(3.10) 


(3.11) 


we  can  reduce  the  condition  (y^,y2)  -  (0,0)  to  (z^,z2)  -  (0,0).  Applying 
Eqs.  (3.10)  and  (3.11)  to  Eq.  (3.9),  the  multiplicative  function  becomes 


F(t) 


-J° 

J0 


Dz(t)  (exp[(im/2K)  I  Dz(t)  [(z1 


f 

J0 


2  2. 

-«  Zl  ) 


+  (z; 


•J  2  2, 

■  3w  z2) 


(3.12) 
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(3.13) 


In  Eq.  (3.12),  J  becomes  unity. 

If  the  action  is  separated  into  the  functionals  with  only  same 
variables  in  the  path  integral,  then  this  integral  can  be  represented  by  the 
multiplication  of  path  integrals  with  each  variable.  Therefore,  Eq.  (3.12) 
becomes 


F(t)  -  FjU)  F2(t) 


8 


-  (J*°  Dz^(t)  exp[(im/2>()  |  dt  (z2-u>2z2) 


x  (J  Dz2(t)  exp((im/2>()  f  dt  (z2-3w2z2)]) 


(3.14) 


Since  F^(t)  and  F2(t)  are  the  path  integrals  of  the  harmonic  oscillator,  the 
evaluation  of  Eq.  (3.14)  gives 


mcj  _ _ 1  h 

2xi)rf  sin(<jt)  sin(,/3a>t)J 


(3.15) 


Hence,  the  propagator  of  DCHO  can  be  written  as 


K(x1,x2,t;  x^.x^.O)  -  [sin(wt) ^sin(/3wt)]  6 


,  iS  /y. 
h  _  c" 


(3.16) 


4.  Propagator  for  the  coupled  harmonic  oscillators 

To  evaluate  the  exact  propagator  expressed  by  Eq.  (3.16),  we  should 
first  obtain  the  propagator  for  CHO.  The  classical  action  of  CHO  is 


a  ,ni/,2  »2  .  2.2  2  .  . 

dr  <2(xc1+xc2)  -  mu  (xcl-xclxc2+xc2))  , 


(4.1) 


where  x^  and  x^  are  the  classical  path  and  velocity,  respectively. 
Integrating  Eq.  (4.1)  over  the  time,  we  get 


Sc  “  2(xclxcl+xc2x 


xc2)lo  *  L 


dr  2  (2xcrxc2)} 
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i: 


««  2 

dr  2  Xc2lXc2^  (2xc2-Xcl)5 


"  2  [xcl(t)Xcl(t)+Xc2(t)Xc2(t)'Xcl(0)Xcl(0)'Xc2(0)Xc2(0) 1  '  (4’2) 

Here  the  second  and  third  terms  become  zero  because  of  the  equations  of 
motion  [see  Eqs.  (2.7)  and  (2.8)],  given  as 


+  u  (2x^-x2)  -  r> 


(4.3) 


x^  +  cj  (2x2’X^)  -  0 


(4.4) 


To  obtain  the  exact  expression  of  Eq.  (4.2),  we  solve  Eqs.  (4.3)  and  (4.4) 
to  obtain 


x^  -  x^(t)  -  A  sin(wt)  +  B  cos(uit)  +  C  sin(/3wt)  +  D  cos(/3u>t)  (4.5) 

X2  ”  X2^C^  ”  A  si-n(wC)  +  B  cos(wt)  -  C  sin(/3u>t)  -  D  cos(/3u>t)  ,  (4.6) 

■  • 

and  x^  and  Xj  are  given,  respectively,  by 

-  x^(t)  -  u)(A  cos(wt)  -  B  sin(wt)  +  /3C  sin(/3wt) 

-  J2D  cin(/3ut))  (47) 


X2  ”  fos(ut)  -  B  sin(wt)  -  73C  cos(/3wt) 


+  /3D  sin(/3wt) ) 


(4.8) 


10 


Equations  (4.5)-(4.8)  give 

x^  -  x^O)  -  B  +  D  (4.9) 

x£  -  x2(0)  -  B  -  D  (4.10) 

x|  -  3^(0)  -  u(A  -i  /3C)  (4.11) 

x'  -  x2(0)  -  cj(A  -  /3C)  (4.12) 

The  time -dependent  constants  A,  B,  D  and  D  obtained  from  Eqs.  (4.5)  and 
(4.6),  and  Eqs.  (4.9)  and  (4.10)  can  be  expressed  as 

A  -  [2  sin(wt) ]  {x^  +  x2  -  (x^  +  x2>  cos(u)t) )  (4.13) 

B  “  2  (xi  +  x2)  (4-14) 

C  -  [2  sin(/3wt ) ]  (x^  -  x 2  +  (x|  -  x£)  cos(/3wt)l  (4.15) 

D  -  2  (x£-x£)  (4.16) 

Substitution  of  Eqs.  (4.5)  -(4.16)  into  (4.2)  gives  the  classical  action: 

Sc  -  {(x2  +  x2  +  x^2  +  Xj2) (cot(wt)  +  /3  cot(/3u>t)]) 


+  2(x^x2  +  x^x2)  (cot(wt)  -  73  cot  (75«t)] 
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-  2(x^x^  +  (  (l/sin(wt)  +  /3/sin(,/3ut)  ]  i 

+  2(x^x£  +  X2xp  [-l/sin(ut)  +  (/3/sin(,/3u>t)  ]  (4.17) 

Combining  Eqs.  (4.17)  and  (3.16),  we  obtain  the  propagator  for  CHO: 

K(x^ ,  x2 ,  t ;  x^.Xj  ,0)  “  [/3/sin(wt)  sin(,/3wt)]‘ 

x  exp(  ( imu)/4h)  [  (x^  +  x^  +  x^  +  x£^) [cot(wt)  +  J2  cot(,/3ojt)  ] 

+  2(x^2  +  xjx^)  [cot(ut)  -  /3  cot(/3a>t)  ]  -  2(x^x^  +  X2x£ ) [ l/sin(wt) 

+  /3/sin(/3wt) ]  +  2(x^x£  +  X2Xp  [-l/sin(u>t)  +  /3/sin(/3iat)  ]  } 

(4.18) 

5 .  Propagator  for  driven  coupled  harmonic  oscillators 

When  we  set  f^(t)  -  f2(t)  -  0,  DCHO  reduces  to  CHO,  whereby  we  can 
write  the  propagator  for  DCHO  as 

K(x1,x2,t;  x^.x^.O)  -  exp[a(t)x^  +  b(t)x1x2  +  c(t)x2  +  d(t)x1 

+  g(t)x2  +  h( t) ]  (5.1) 

Here  a(t),  b(t) ,  c(t),  d(t) ,  f(t)  and  h(t)  are  time -dependent  functions 
including  xj^  and  x£ ,  which  need  to  be  determined.  Equation  (5.1)  must 
satisfy  the  Schroedirger  equation 
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iJrt(3K/3t)  -  H  K  (5.2) 

Substitution  of  Eq.  (5.1)  into  Eq.  (5.2)  gives  the  time -dependent 
coefficients 


a(t)  -  jjj  [4a2 (t)  +  c2(t)]  +  mw2/i *  (5.3) 
b(t)  -  [4b2 (t)  +  c2(t)}  +  mw2/itf  (5.4) 
c(t)  -  ^  [a(t)c(t)  +  b(t)c(t)]  -  mcj2/iy  (5.5) 
d(t)  -  [2a(t)d(t)  +  C(t)g(t) ]  +  (j)fx(t)  (5.6) 
g(t)  -  ^  [2b(t)g(t)  +  c(t)d(t) ]  +  (J)f2(t)  (5.7) 
h(t)  -  &  [d2(t)  +  g2(t)  +  2a(t)  +  2b(t) ]  .  (5.8) 


Since  Eqs .  (5.3)  and  (5.4)  have  the  same  form,  we  get 

a(t)  -  b(t)  .  (5.9) 

Substituting  Eq.  (5.9)  into  Eq.  (5.5)  and  changing  the  variables  a  and  c 
into 


r)  —  a  +  c/2  , 

(5.10) 

f  -  a  -  c/2 

(5.11) 

1 


we  obtain  two  ordinary  differential  equations: 


•  2iH  2  _  ma> 

*  “  m  *»  +  2i* 


(5.12) 


:  m  r2  3mo £ 

f  m  f  2i# 


(5.13) 


The  solutions  of  Eqs .  (5.12)  and  (5.13)  are  given  by 


V  ~  ~ ^  cot(«t  +  9 


(5.14) 


f  -  ^2^"m  c°c(y3<*>t  +  «2) 


(5.15) 


where  9^  and  #2  are  the  constants  to  be  determined.  The  time -dependent 
coefficients  a(t) ,  b(t)  and  c(t)  obtained  in  comparison  with  Eqs.  (5.10), 
(5.11),  (5.14)  and  (5.15)  are  given  as 


i(t)  -  b(t)  -  [cot(wt  +  9:)  +  J1  cot(/3wt  +  9n) 

4  ft  i  Z 


(5.16) 


c(t)  -  [cot(wt  +  9^)  -  73  cot(/3wt  +  9 y) 


(5.17) 


Equations  (5.16)  and  (5.17)  do  not  include  the  driven  forces  f^(t)  and 

f 2 ( t ) .  Therefore,  through  setting  f^(t)  -  f2(t)  -  0,  Eqs.  (5.16)  and  (5.17) 

2  2 

do  not  change  at  all  and  should  be  equal  to  the  coefficients  of  x^  and  x2  in 
Eq.  (4.18).  Comparison  of  these  two  equations  shows  9 ^  and  §2  to  be  zero. 
Substituting  Eq.  (5.9)  into  Eqs.  (5.6)  and  (5.7)  and  changing  variables  d 
and  g  into 


14 


rt 

p  -  [  l/sin(ut)  ]  {J  dr  ^  [f^r)  +  f  2  ( »- )  ]  sin(u>t)  +  a}  (5.22) 

ft 

a  -  [ l/sin(/3wt) ]  (j  dr  J  [f1(r>  -  f2<r)]  sin(73wt)  +  0)  (5.23) 

where  a  and  0  are  constants  to  be  determined.  We  can  obtain  the  time- 
dependent  coefficients  d(t)  and  g(t)  by  substituting  Eqs.  (5.22)  and  (5.23) 
into  Eqs.  (5.18)  and  (5.19): 

d(t)  -  [i/2)(sin(wt) ]  J  dr(f^(r)  +  f2(r)]  sin(wr) 

+  (i/2)(sin(y3«t)]  J  dr[f1(r)  -  f 2 ( r )  ]  sin(73o>r) 

+  [a/2sin(wt)]  +  [£/2sin(./3ut)  ] 


(5.24) 
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g(t)  -  ^i/2Xsin(wt)j  |  dr  [f^r)  +  f 2 ( r )  ]  sin(u>r) 

-  j^i/2Xsin(73wt)j  J  dr  (f^r)  -  €2 ( r )  ]  sin(/3ur) 


+  [a/2sin(wt) ]  -  [£/2sin(-/3wt)  ] 


(5.25) 


Substitution  of  Eqs.  (5.16),  (5.17).  (5.24)  and  (5.25)  into  Eq.  (5.8)  yields 


h(t)  -  -  (o2cot(cjt)  +  (02//3)  cot(/3wt)] 


[a/musLn(uC) ] 


lo  1£l 


(r)  +  fj(r  )  ]  sin[a)(t-r  )  ] 


ft 

-  [0//3mwsin(/3cjt)  1  J  dr  (f^r)  -  f2(r)]  sin[/3w(t-r; 


[l/4iXmwsin(wt)  ]  [  dr  f  di/  [  f..  ( r )  +  f  9  ( r )  ] 
J0  J0 


x  [^(v)  +  f2(*')l  sin[«(t-r)]sin(wi/) 


[l/4,/3iXmw  sin(/3wt)]  [  dr  [  d»/  [f,(r)  -  f_(r)] 

J0  J0 


x  ( f  1(ix)  -  f2(l/)]sin[73uj(t-r)]  sin(/3w»/)  -  in[sin(ut) 


x  sin(/3wt)  +  S 


(5.26) 


Here,  S  is  also  a  constant  to  be  determined.  When  setting  f^(t)  -  f2(t) 
0,  Eqs.  (5.24)  and  (5.25)  should  be  reduced  to  the  coefficients  of  and 
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,  and  Eq.  (5.26)  should  also  be  reduced  to  the  terms  in  the  exponent  in 
Eq.  (4.18).  Comparison  between  them  gives  the  constants  a,  £>  and  &: 


a  “  I>(  (xl  +  x2^ 


(5.27) 


“  (xi  -  x2>  ■ 


(5.28) 


-k 

,  .  ,  3  ma). 

5  ~  in(2iriK) 


(5.29) 


Substitution  of  the  above  results  into  Eq.  (5.1)  gives  the  propagator  for 


DCHO : 


K(xlfx2,t;  x£,x£,0)  -  ^j-jj|,/3/(sin(wt)  sin(/3wt)  ] 
x  exp^®jj|  |^(x2  +  x2  +  x^2  +  x2^)[cot(wt)  +  j3cot(,JZiot) 


+  2(x^x2  +  x|Xj) (cot(ut)  -  73cot(/3wt)]  -  2(x^x^  +  x^x^) ( l/sin(wt) 


+  j3/stn(/3at) 


+  2(x^x2  +  x^x2) [ -l/sin(ut)  +  /3/s in (/Jut ) ] 


+  ( ( l/sin(wt)  ]  J  dr  (f^r)  +  f  2  (r )  ]  :,in(u>r) 


[l/sin(/3»t)]  J  dr  |fx<r)  -  f 2 ( r ) ]  sin(/3wr)} 
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+  m T  Ul/*in(»t)] 


I  dr  ffl 
J0 


(r)  +  f j ( r ) ]  sin (ur) 


-  I  L/sin(/3wt)  ]  |  dr  (f^Cr)  -  f2(r>]  sin(/3wr>) 

+  ^  ( [l/sin(wt) ]  J  dr  ( f ^ ( r )  +  f 2 ( r ) ]  sin[w(t-r)] 

+  [l/sin(/3ut) ]  J  dr  [ f ^ ( r )  -  C  ^ ) ]  sin(/3wr)) 


Av  * 


f[l/sin(wt)]  J  dr  [f^(r)  +  f 2 ( r ) ]  sin[w(t-r)] 


-  [l/sin(/3u>t)  ]  J  dr  (f^r)  -  f2(r))  sin(/3ur)) 


,22  rc  rt 

[1/m  <j  sin(wt)  ]  dr  di/  [f,(r)  +  f 0 ( r )  ]  [ f..  (is)  +  f„ 

Jq  Jq  1  1 


(*>] 


x  sin[w(t-r)]  sin(wi/) 


-  [l/y3m2a;2sin(/3wt)  J  J  dr  J  di/  [ f x ( r )  -  f2(r)  ]  -  f 2 C»^ )  I 


x  sin[/3w(t-r )  ]  sin(/3w»/) )  j  J  . 


C 


t 


.30) 
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6 .  Energy  expectation  values  of  coupled  hanonic  oscillators 

The  Hamiltonian  of  CHO  is 


H  “  2m  (P1  +  p2}  +  “u,2^xi  '  xix2  +  x2} 


(6.1) 


Using  Eqs.  (3.1)  and  (3.2)  with  Eq.  (6.1),  we  obtain  the  Schroedinger 
equation. 


iHO/3t)  0(x1,x2,t)  -  H  \i(x1,x2,t) 


(6.2) 


where  H  is  the  Hamiltonian  operator  in  which  the  momentum  is  changed 
into  pt  -  (H/i)(3/3xi) .  Since  Eq.  (6.2)  can  be  separated  into  time  and 
coordinate  parts,  we  may  write 


-iH  t/>( 
K(t)  -  e  op 


(6.3) 


Hop'iin>  "  Ein^>n>  •  (i,n  -  1,2,3,...) 


(6.4) 


Here  the  states  |i,n>  are  the  complete  set  with  energy  eigenvalues  of  H  . 
Since  the  function  with  states  |2,n>  can  be  expressed  by 


^in^Xl,X2)  -  <xi*x2li*n> 


(6.5) 


the  propagator  at  t  >  0  becomes 


K(x^,x2,t;  x^.x^.O)  -  <x1,x2|e  p  | 


iH  t/)i 

op  |x',x'> 
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r-  r-  r-  r-  -  iH 

~)  )  )  /  <x1,x2|in>  <in|e  0p  | i' ,n'>  <i’  ,n' |x|,x£> 

i  n  V  n' 

-}  2  *in(xrx2>  e  ^n<xrx2>  ' 

l  n 


Equation  (6.6)  should  be  the  same  as  Eq.  (4.18).  Setting  and  x£  - 

X£  in  Eq.  (4.18)  and  integrating  over  x^  and  x2 ,  we  get 


2  2  J  j  dXldX2  ^n(xl-x2)  6  V(xl'x2) 

l  n 


(6.7) 


and 

{  J  dxldx2  ^{^/[  sin(wt)sin(y3u>t) 

x  exp(-^p[  (x^  +  x2)2  -  /3(x^  -  Xj)2][cot(wt)  -  l/sin(wt)  ]  ) 

-  -  ^  [ sin(wt/2)sin(73/2wt) ]  ^  (6.8) 

Hence ,  we  have 

r*  r"  iE »  t/)£  -  _  1 

Me  n  -  -  j  [sin(wt/2)sin(73/2wt) ] 

Z  n 

-  (e-^/d  -  ea„t)H,.iy5,t/2/a  .  .-t/s*),. 
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-  ^  ^  exp{  -  iu)t[  (2  +  +  /3(n  +  2)])  (6.9) 

2-0  n-0 

Therefore  the  expectation  values  of  CHO  becomes 

E2n  “  [(i  +  b  +  ^(n  +  2)]  ■  (6-10) 

7 .  Results  and  discussion 

In  the  previous  sections  we  have  obtained  the  exact  propagators  [Eqs. 

(4.18)  and  (5.30)]  for  CHO  and  DCHO  by  the  path- integral  method.  The  forms 
of  the  propagators  are  new.  Setting  f(t)  -  0,  Eq.  (5.30)  is  reduced  to  Eq. 

(4.18) .  Although  DCHO  is  a  nonconservative  system,  the  quantum -mechanical 

problem  for  the  momentum  operator  does  not  appear  because  the  canonical 

13 

momentum  is  equal  to  the  kinetic  momentum  in  our  derivation. 

Making  use  of  Eq.  (4.18),  we  have  obtained  the  energy  expectation 

values  [Eq.  (6.10)]  for  CHO,  given  by  the  sum  of  two  energy  expectation 

values  corresponding  to  the  quantum  states  of  two  oscillators.  Even  though 

we  have  not  evaluated  the  wavefunction  of  CHO,  we  may  easily  surmise  that 

the  wavefunction  will  be  given  by  the  multiplication  of  two  wavefunctions 

for  two  oscillators.  In  the  case  of  DCHO,  one  cannot  easily  apply  Eq. 

(5.20)  to  obtain  the  energy  expectation  values,  since  this  equation  cannot 

be  expressed  in  the  form  of  Eq.  (6.6),  and  one  should  recognize  that  the 

energy  operator  is  not  equal  to  the  Hamiltonian  operator  in  a 

9 

nonconservative  system. 

The  evaluations  for  the  wavefunctions,  energy  expectation  values  for 
CHO  and  DCHO,  and  propagator  and  other  physical  quantities  for  n  coupled  and 
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n  driven  coupled  harmonic  oscillators  (arbitrary  n)  are  in  progress  and  will 
be  reported  in  the  near  future . 


Acknowledgments 

This  research  was  supported  in  part  by  a  grant  to  Korea  University 
from  the  BSRI  Program,  Ministry  of  Education  1987,  Republic  of  Korea,  and  in 
part  by  the  Office  of  Naval  Research,  the  Air  Force  Office  of  Scientific 
Research  (AFSC) ,  under  Contract  F49620-86-C-0009 ,  and  the  National  Science 
Foundation  under  Grant  CHE- 8620274.  The  United  States  Government  is 
authorized  to  reproduce  and  distribute  reprints  for  governmental  purposes 
notwithstanding  any  copyright  notation  hereon. 


22 


References 

1.  R.  P.  Feynman  and  A.  R.  Hibbs ,  "Quantum  Mechanics  and  Path  Integrals" 
(McGraw-Hill,  New  York,  1965). 

2.  J.  T.  Marshall  and  J .  L.  Pell,  J.  Math.  Phys .  22.  1297  (1979);  G.  E. 
Prince  and  C.  J.  Eliezer,  J.  Phys.  A  H,  815  (1980). 

3.  L.  F.  Landovitz,  A.  M.  Levine  and  W.  M.  Schreiber,  Phys.  Rev.  A  20 . 

1162  (1979);  J.  Math.  Phys.  21,  2159  (1980);  A.  M.  Levine,  E.  Ozizmir 
and  W.  M.  Schreiber,  J.  Chem.  Phys.  2&.  291  (1983). 

4.  A.  C.  Rose-Innes  and  E.  H.  Rhoderick,  "Introduction  to 
Superconductivity"  (Pergamon,  Oxford,  1978). 

5.  V.  B.  Braginsky  and  Yu.  I.  Vorontsor,  Usp.  Fiz,  Nauk  114 .  41  (1974);  U. 
G.  Unruh,  Phys.  Rev.  D  M,  1764  (1978);  12,  2888  (1979). 

6.  R.  King,  B.  A.  Huberman  and  J.  Barchas ,  in  "Synergetics  of  the  Brain," 
ed.  by  H.  Haken  (Springer,  Berlin,  1983). 

7.  P.  Caldirola,  Nuovo  Cimento  18,  393  (1941);  B  77,  241  (1983);  E.  Kanai, 
Prog.  Theor.  Phys.  3,  440  (1948). 

8.  H.  Dekker,  Phys.  Rep.  80,  1  (1981);  B.  K.  Cheng,  J.  Phys.  A:  Math.  Gen. 
12.,  2475  (1984);  K.  H.  Yeon,  C.  I .  Um  and  T.  F.  George,  Phys.  Rev.  A 
36,  5287  (1987);  E.  Eckhardt,  Phys.  Rev.  A  35,  5191  (1987). 

9.  V.  V.  Dodonov  and  V.  I.  Manko,  Phys.  Rev.  A  20,  550  (1979). 

10.  I.  R.  Senitzky,  Phys.  Rev.  H£,  670  (1960). 

11.  D.  M.  Greenberger,  J.  Math.  Phys.  22,  762  (1979). 

12.  J.  J.  Cervero  and  J.  Villaroel,  J.  Phys.  A:  Math.  Gen.  12,  2963  (1984). 

13.  C.  I.  Um,  K.  H.  Yeon  and  W.  H.  Kahng,  J.  Phys.  A:  Math.  Gen.  22,  611 

(1987);  J.  Korean  Phys.  Soc.  12,  1,  8  (1986);  New  Phys.  (Korean  Phys. 
Soc.)  22,  100,  117  (1985). 


CL/1113/86/2 


TECHNICAL  REPORT  DISTRIBUTION  LIST.  GEN 


No. 

Copies 

Office  of  Naval  Research  2 

Attn:  Code  1113 
800  N.  Quincy  Street 
Arlington,  Virginia  22217-5000 

Or.  8ernard  Ocuda  1 

Naval  Weapons  Support  Center 
Code  50C 

Crane,  Indiana  47522-5050 


Naval  Civil  Engineering  Laboratory  1 
Attn:  Or.  R.  W.  Drisko,  Code  L52 
Port  Hueneme,  California  93401 


Defense  Technical  Information  Center  12 


Building  5,  Cameron  Station  high 

Alexandria,  Virginia  22314  quality 

OTNSRDC  1 

Attn:  Or.  H.  Slngerman 
Applied  Chemistry  Division 
Annapolis,  Maryland  21401 

Or.  William  Tolies  1 

Superintendent 


Chemistry  Division,  Code  6100 
Naval  Research  Laboratory 
Washington,  D.C.  20375-5000 


No. 

Copies 

Or.  David  Young  1 

Code  334 

NOROA 

NSTl,  Mississippi  39529 

Naval  Weapons  Ceiter  1 

Attn:  Dr.  Ron  Atkins 
Chemistry  Olvlslon 
China  Lake,  California  93555 

Scientific  Advisor  1 

Commandant  of  the  Marine  Corps 
Code  RD-l 

Washington,  D.C.  20380 

U.S.  Army  Research  Office  ] 

Attn:  CRD-AA-IP 
P.0.  Box  12211 

Research  Triangle  Park,  NC  27709 

Mr.  John  Boyle  1 

Materials  8ranch 
Naval  Ship  Engineering  Center 
Philadelphia,  Pennsylvania  19112 

Naval  Ocean  Systems  Carter  1 

Attn:  Dr.  S.  Yamamoto 
Marine  Sciences  Division 
San  Diego,  California  91232 

Dr.  David  L.  Nelson  1 

Chemistry  Division 
Office  of  Naval  Research 
800  North  Quincy  Street 
Arlington,  Virginia  22217 


1 


OL/1113/86/2 


ABSTRACTS  DISTRIBUTION  LIST,  056/625/629 


Dr.  J.  E.  Jensen 
Hugnes  Research  Laboratory 
3011  Malibu  Canyon  Road 
Malibu,  California  90265 

Dr.  J.  H.  Weaver 

Department  of  Chemical  Engineering 
and  Materials  Science 
University  of  Minnesota 
Minneapolis,  Minnesota  55455 

Dr.  A.  Reisman 

Microelectronics  Ceiter  of  North  Carol 
Research  Triangle  Park,  North  Carolina 
27709 

Dr.  M.  Grunze 

Laboratory  for  Surface  Science  and 
Technology 
University  of  Maine 
Orono,  Maine  04469 

Dr.  J.  Butler 

Naval  Research  Laboratory 

Code  6115 

Washington  D.C.  20375-5000 

Dr.  L.  Interante 
Chemistry  Department 
Rensselaer  Polytechnic  Institute 
T  roy ,  New  York  12181 

Dr.  Irvin  Heard 

Chemistry  and  Physics  Department 
Lincoln  University 

Lincoln  University,  Pennsylvania  19352 

Dr.  K.J.  Klaubunde 
Department  of  Chemistry 
Kansas  State  University 
Manhattan,  Kansas  66506 


Dr.  C.  8.  Harris 
Department  of  Chemistry 
University  of  California 
Berkeley,  California  94720 


Dr.  F.  Kutzler 
Department  of  Chemistry 
Box  5055 

Tennessee  Technological  University 
Cookesville,  Tainessee  38501 

na  Dr.  D.  Di Leila 

Chemistry  Department 
George  Washington  University 
Washington  D.C.  20052 

Or.  R.  Reeves 
Chemistry  Department 
Renssaeler  Polytechnic  Institute 
Troy,  New  York  12181 

Dr.  Steven  M.  George 
Stanford  University 
Department  of  Chemistry 
Stanford,  CA  94305 

Dr.  Mark  Johnson 
Yale  University 
Department  of  Chemistry 
New  Haven,  CT  06511-8118 

Or.  W.  Knauer 

Hughes  Research  Laboratory 
3011  Malibu  Canyon  Road 
Malibu,  California  90265 


6 


a/1113/86/2 


ABSTRACTS  DISTRIBUTION  LIST,  056/625/629 


Cr.  G.  A.  Somorjal 
Department  of  Chemistry 
University  of  California 
Berkeley,  California  94720 

Or.  J.  Murday 

Naval  Research  Laboratory 

Code  6170 

Washington,  D.C.  20  3  75  -  5000 

Dr.  J.  B.  Hudson 
Materials  0  Vision 
Rensselaer  Polytechnic  Institute 
Troy,  New  York  12181 

Or.  Theodore  E.  Madey 
Surface  Chemistry  Section 
Department  of  Commerce 
National  Bureau  of  Standards 
Washington,  D.C.  20234 

Dr.  J.  E.  Oemuch 
IBM  Corporation 

Thomas  J.  Watson  Research  Center 
P.O.  Box  218 

Yorktown  Heights,  Nor  York  10598 

Dr.  M.  G.  Lagally 
Department  of  Metallurgical 
and  Mining  Engineering 
University  of  Wisconsin 
Madison,  Wisconsin  53706 

Dr.  R.  P.  Van  Duyne 
Chemistry  Department 
Northwestern  University 
Evanston,  Illinois  60637 

Or.  J.  M.  White 
Department  of  Chemistry 
University  of  Texas 
Austin,  Texas  78712 

Dr.  O.  E.  Harrison 
Department  of  Physics 
Naval  Postgraduate  School 
Monterey,  California  93940 


Dr,  R.  L.  Park 

Director,  Center  of  Materials 
Research 

University  of  Maryland 
College  Park,  Maryland  20742 

Dr.  W.  T.  Perla 

Electrical  Engineering  Department 
University  of  Minnesota 
Minneapolis,  Minnesota  55455 

Dr.  Keith  H.  Johnson 
Department  of  Metallurgy  and 
Materials  Science 

Massachusetts  Institute  of  Technology 
Cambridge,  Massachusetts  02139 

Dr.  S.  Sibener 
Tepartment  of  Chemistry 
v.  .ies  Franck  Institute 
5640  Ellis  Avenue 
Chicago,  Illinois  60637 

Dr.  Arnold  Green 

Quantum  Surface  Dynamics  Branch 

Code  3817 

Naval  Weapons  Canter 
China  Lake,  California  93555 

Dr.  A.  Wold 

Department  of  Chemistry 
Brown  University 
Providence,  Rhode  Island  02912 

Dr.  S.  L.  Bernasek 
Department  of  Chemistry 
Princeton  University 
Princeton,  New  Jersey  08544 

Dr.  W.  Kohn 
Department  of  Physics 
University  of  California,  San  Diego 
La  Jolla,  California  92037 


7 


01/1113/86/2 


ABSTRACTS  DISTRIBUTION  LIST,  056/625/629 


Dr.  F.  Carter 
Code  6170 

Naval  Research  Laboratory 
Washington,  O.C.  20375-5000 

Or.  Richard  Colton 
Code  6170 

Naval  Research  laboratory 
Washington,  D.C.  20375-5000 

Dr.  Dan  Pierce 
National  Bureau  of  Standards 
Optical  Physics  Oivision 
Washington,  D.C.  20234 


Dr.  R.  Stanley  Williams 
Department  of  Chemistry 
University  of  California 
Los  Angeles,  California  90024 

Dr.  R.  P.  Messmer 
Materials  Characterization  Lab. 
General  Electric  Company 
Schenectady,  New  York  22217 

Dr.  Robert  Gomer 
Department  of  Chemistry 
James  Franck  Institute 
5640  Ellis  Avenue 
Chicago,  Illinois  60637 

Dr.  Ronald  Lee 
R301 

Naval  Surface  Weapons  Center 
White  Oak 

Silver  Spring,  Maryland  20910 

Dr.  Paul  Schoen 
Code  6190 

Naval  Research  Laboratory 
Washington,  D.C.  20375-5000 


Dr.  John  T.  Yates 
Department  of  Chemistry 
University  of  Pittsburgh 
Pittsburgh,  Pennsylvania  15260 

Dr.  Richard  Greene 
Code  5230 

Naval  Research  Laboratory 
Washington,  O.C.  20375-5000 

Dr.  L.  Kesmodel 
Department  of  Physics 
Indiana- Uni versity 
Bloomington,  Indiana  47403 

Dr.  K.  C.  Janda 
University  of  Pittsburg 
Chemistry  Building 
Pittsburg,  PA  15260 

Dr.  E.  A.  Irene 
Department  of  Chemistry 
University  of  North  Carolina 
Chapel  Hill,  North  Carolina  27514 

Dr.  Adam  Heller 

Bell  Laboratories 

Murray  Hill,  Nav  Jersey  07974 

Or.  Martin  Fleischmann 
Department  of  Chemistry 
University  of  Southampton 
Southampton  509  5NH 
UNITED  KINGDOM 

Dr.  H.  Tachikawa 
Chemistry  Department 
Jackson  State  University 
Jackson,  Mississippi  39217 

Dr.  John  W.  Wilkins 
Cornell  University 
Laboratory  of  Atomic  and 
Solid  State  Physics 
Ithaca,  Nav  York  14853 


8 


01/1113/86/ 2 


ABSTRACTS  DISTRIBUTION  LIST,  056/625/629 


Or.  R.  G.  Wallis 
Department  of  Physics 
University  of  California 
Irvine,  California  92664 

Or.  0.  Ramaker 
Chemistry  Department 
George  Washington  University 
Washington,  O.C.  20052 

Or.  J.  C.  Hemminger 
Chemistry  Department 
University  of  California 
Irvine,  California  92717 

Dr.  T.  F.  George 
Chemistry  Department 
University  afJ*t5chester 
RochestarV^ew  York  14627 

1h\  G.  Rubloff 
IBM 

Thomas  J.  Watson  Research  Center 
P.0.  Box  218 

Yorktown  Heights,  N»*  York  10598 

Dr.  Horia  Metiu 
Chemistry  Department 
University  of  California 
Santa  Barbara,  California  93106 

Dr.  W.  Goddard 

Department  of  Chemistry  and  Chemical 
Engineering 

California  Institute  of  Technology 
Pasadena,  California  91125 

Dr.  P.  Hansma 
Department  of  Physics 
University  of  California 
Santa  Barbara,  California  93106 

Dr.  J.  Baldeschwieler 
Department  of  Chemistry  and 
Chemical  Engineering 
California  Institute  of  Technology 
Pasadena,  California  91125 


Dr.  J.  T.  Keiser 
Department  of  Chemistry 
University  of  Richmond 
Richmond,  Virginia  23173 

Dr.  R.  W.  Plummer 
Department  of  Physics 
University  of  Pennsylvania 
Philadelphia,  Pennsylvania  19104 

Dr.  E.  Yeager 
Department  of  Chemistry 
Case  Western  Reserve  University 
Cleveland,  Ohio  41106 

Dr.  N.  Winograd 
Department  of  Chemistry 
Pennsylvania  State  University 
University  Park,  Pennsylvania  16802 

Dr.  Roald  Hoffmann 
Department  of  Chemistry 
Cornell  University 
Ithaca,  New  York  14853 

Dr.  A.  Steckl 

Department  of  Electrical  and 
Systems  Engineering 
Rensselaer  Polytechnic  Institute 
Troy,  N ewYork  12181 

Dr.  G.H.  Morrison 
Department  of  Chemistry 
Cornell  University 
Ithaca,  New  York  14853 


